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Abstract 

The role of magnetic shear for zonal flow generation by ion-temperature- 
gradient (ITG-) and trapped electron (TE-) mode turbulence is stud- 
ied analytically using fluid descriptions. The scaling of the zonal flow 
(ZF) growth rate with magnetic shear is examined and compared with 
linear growth rates for typical tokamak parameter values. The results 
indicate that large levels of ZF are obtained in regions of negative mag- 
netic shear, in particular for ZF driven by TE mode turbulence. The 
strong magnetic shear scaling obtained for TE mode driven zonal flows 
originates from the bounce average of the electron magnetic drifts. 
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1 Introduction 



The study of turbulence suppression and generation in magnetically confined 
fusion plasma research has attracted much interest recently. It is well estab- 
lished that the interaction between small scale turbulence and zonal flows, 
which are radially localized flows propagating mainly in the poloidal direc- 
tion, can reduce the radial transport by shearing the eddies of the driving 
background turbulence pQ- [3] . Accordingly, the self-consistent generation of 
zonal flows by non-linear interactions among drift waves has recently been 
studied both analytically [5]- [H] and numerically [15]- |23j. 

The purpose of the present letter is extend previous analytical studies [H]- 
[TU] to investigate the effects of magnetic shear on the zonal flow generation 
by ion-temperature-gradient (ITG-) and trapped electron (TE-) mode turbu- 
lence, respectively. The ITG- and TE-mode driven cases will be considered 
separately and the corresponding ZF generation by ITG and TE mode tur- 
bulence will be calculated and compared. 

The fluid models used in the present Letter have been extensively de- 
scribed in previous papers and are only briefly summarized here. For the 
ITG mode description, a reduced version of the standard Weiland model [23] 
is utilized. In this formulation the ITG model includes the ion continuity 
and the ion temperature equations, neglecting the influence of parallel ion 
momentum and electromagnetic effects, while the electrons are assumed to 
be Boltzmann distributed. The TEM model is based on the bounce averaged 
electron continuity and electron temperature equations for the trapped elec- 
tron fluid [21]. The ion and trapped electron mode fluid models are fairly 
symmetric except for the FLR terms which are included in the ion fluid model 
and the bounce average which is performed on the electron quantities. The 
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equations are: 
dn 



dt 



+ V ■ (njV E + njiJ+j) + V ■ (njVpj + njV nj ) = (1) 
3 dT 

+ • tT,. + V ■ <g = (2) 

where rij, Tj are the density and temperature perturbations (j = i and j = et 
represents ions and trapped electrons) and Vj = ve + v* + vpj + v n j, ve is the 
E x B velocity, v+ is the diamagnetic drift velocity, vpj is the polarization 
drift velocity, v n j is the stress tensor drift velocity and qj is the heat flux. The 
derivative is defined as d/dt = d/dt + p s c s zx V0- V and <fi is the electrostatic 
potential. The linear solutions to Eqs 1 - 3 for TE modes and ITG modes 
respectively are given by 

UrTEM = -y (t (1 - e n g e ) - Y € n9e^j (4) 

Item = k y ^e n g e (r] e - r} eth ) (5) 
2 £ 10 , ie n g e £ 
Veth = 3-2 + 9f^ + ^ + 4^ (6) 



iitg = 1 ^ k 2 xl^r-to - vuh) (8) 



2 T T (T 10T\ 

* 3-2 + 4^- + e ^U + — J (9) 
where uj = u r + ry. Here, e n = L n /R, rjj = L n /L T j, r = T e /Tj, ctj = 
^ (1 + ?7i), £ is the fraction of trapped electrons and gj (j = i, e) are defined 
below. 

The zonal flow evolution is treated by the vorticity equation. The gener- 
ation of ZF by ITG- and TE-modes are studied separately for rji >> i] e (ITG 
dominated case) and r\ e >> rji (TEM dominated case) respectively. In these 
limits it is assumed that the electron and ion responses are decoupled and 



hence the effects of ion perturbations on the TE mode are neglected (and 
vice versa). The previous work on zonal flow generation by ITG and TE 
modes [9]- [10] is here extended to include the magnetic shear dependence 
through the magnetic drift frequencies and the perpendicular wave vector. 
The effect of magnetic shear enters through the averaged magnetic drift fre- 
quency for the electrons and ions. The ion fluid model is averaged using 
an approximate eigen-mode function resulting in an averaged magnetic drift 
frequency: ip — 7^(1 + cos(0)) and \9\ < n, giving oc e n (| + |s — e) and 
(k\) oc (1 — 2.5s 2 + ^s 2 ) [25J. Here s is the magnetic shear and e = a/R 
is the inverse aspect ratio parameter. The effects on zonal flow generation 
of magnetic shear have been treated for ITG modes in Ref. [T3] using the 
reductive perturbation method for the zonal flow generation. 

The electron fluid equations are bounce averaged and the electron mag- 
netic drift is replaced by the precession frequency of trapped electrons (uJDe) = 
^Deo9e where uJDeo = 2ke ^ sCs and g e — \ + |s [26J. The original formulation 
of the Weiland model for deeply trapped electrons is recovered for s = 9/8 
in this model. 

The zonal flow dispersion relations are found using the vorticity equation 
and combining this with the wave kinetic equation as an additional relation 
between the background turbulence and the zonal flow. The average of the 
ion continuity equation over the magnetic surface and over fast small scales 
employing the quasi-neutrality the evolution of the mean flow is obtained 

- ^ ={ l+s ?) v * + ^ (=*£*•) 

where it is assumed that only the small scale self interactions are the im- 
portant interactions in the RHS. Here 5i m is the Kronecker delta with m = 
{i,et}. It is also assumed that there is a sufficient spectral gap between 
zonal flows and the background turbulence. In the case of zonal flows (fl) 



generated from ITG mode turbulence the dispersion relation is [H] 

(n + tfiql) (Q - q x v gx f = -ql(l+ l - + U^j |0 fc | 2 fi (11) 
and in the TEM generation case we obtain [10] 



Q = iq x k y \J \<p k \ 2 . (12) 

Here, Q is the zonal flow growth rate, k x , k y is the wave numbers for the 
background turbulence, q x is the wavenumber of the zonal flow, v gx is the 
group velocity for the drift wave and /j, is a corresponding collisional damping 
(which is neglected in this study). The 5 is a wave number independent factor 
given by 

A k ky ( 2 /_ 1 

K / i 



5 = Aft%^-3( 1 + -) e ^) (13) 



A fc = '-j- \ l-e ngi +^-e n gA . (14) 

Hence, the zonal flow growth rate scales as Q oc \<pk\- Proceeding from 
here, several different assumptions are possible on the saturation level for 
the background turbulence in absence of zonal flows. In the present work 
only the mixing length saturation level (Tj = k Lt ) is considered. 



2 Results 

The dispersion relations (Eqs 1 and 2) for the zonal flow (Q) is solved nu- 
merically. Special attention is given to the magnetic shear dependency of the 
zonal flows generated from ITG mode and TEM turbulence, respectively. 

In Figure 1, the magnetic shear scaling of the zonal flow growth rate driven 
by ITG modes (normalized to the linear ITG growth rate) is displayed. The 
parameters are k x = k y = 0.3, q x = 0.3, r = 1.0 and e = 0. In Figure la 
the other parameters are rj e = (suppressing the TE mode), r\i = 4.0, ft = 
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Figure 1: The ZF growth (normalized to the ITG mode growth) as a function 
of magnetic shear s with e n as a parameter. The results are shown for e n = 1.5 
(rings), e n = 1.0 (asterisk) and e n = 0.5 (diamonds). The other parameters 
are k x = k y = 0.3, q x = 0.3, e = 0, i] e = 0, i]i = 4.0 and fraction of trapped 
electrons f t = 0. 

(fraction of trapped electrons ) and e n is given by e n = 1.5 (rings), e n = 1.0 
(asterisk) and e n = 0.5 (diamond). 

The corresponding result for the TEM driven case is shown in Figure 
2 for f]i = (suppressing the ITG mode), i] e = 4.0 and f t = 0.5. The 
parameter e n is given by e n = 2.0 (plus), e n = 1.5 (rings) and e n = 1.0 
(asterisk). For sufficiently large negative shear the linear TEM is stabilized 
due to the change of sign of the precession frequency. For the zonal flows 
generated by TEM background turbulence there is a significant increase in 
growth rate in the weak and negative shear regions where •jzf >> Item- For 
zonal flows generated by ITG turbulence on the other hand, only a modest 
increase in the corresponding region is found. For positive magnetic shear, 
the ZF scaling with magnetic shear is rather weak. In general we obtain 
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Figure 2: The ZF growth (normalized to the TEM growth) as a function of 
magnetic shear s with e n as a parameter. The results are displayed for e n = 
2.0 (plus), e n = 1.5 (rings) and e n = 1.0 (asterisk). The other parameters 
are k x = k y = 0.3, q x = 0.3, e = 0, r/j = 0, r\ e = 4.0 and fraction of trapped 
electrons f t = 0.5. 

^zf(ITGM) > -f ZF {TEM) in this region. We have verified that the ZF 
growth scales weakly with the other parameters. 

Next, the scaling of zonal flow generation with magnetic shear with nor- 
malized temperature gradient as a parameter is displayed. The parameters 
are k x = k y = 0.3, q x = 0.3, r = 1.0 and e = 0. In Figure 3 the other 
parameters are r\ e = 0, e n = 1.0, ft = with rji = 6.0 (squares) and r\i = 4.0 
(asterisk). The corresponding result for the TEM driven case is displayed in 
Figure 4 for rji = 0, e n = 1.0 and f t = 0.5. The parameter r/ e is given by 
r\ e = 6.0 (squares) and r\ e = 4.0 (asterisk). As found in Figures la and lb for 
weak/negative magnetic shear, a substantial increase in the zonal flow gener- 
ation for TEM turbulence whereas a modest increase is found for ZF driven 
by ITG turbulence. The dependence on normalised temperature gradient is 
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Figure 3: The ZF growth (normalized to the ITG mode growth) as a function 
of magnetic shear s with normalized temperature gradient as a parameter. 
The results are shown for rji = 6.0 (squares) and rji = 4.0 (asterisk). The 
other parameters are k x = k y = 0.3, q x = 0.3, e = 0, rj e = 0, e n = 1.0 and 
fraction of trapped electrons f t = 0. 

weak except close to the linear threshold. 

3 Summary 

In summary, the results indicate that the effects of magnetic shear for ZF gen- 
eration are significantly different for ITG- and TE-mode driven cases. This 
difference is related to the different shear dependence of the ion and electron 
magnetic drift frequency respectively. Hence, for the TE mode dominated 
case, a weak/negative magnetic shear results in a large zonal flow growth rate 
(Izf > Item)- This may result in a strong stabilization of the turbulence. 
For the ITG dominated case on the other hand, the results indicate that the 
effects of magnetic shear on the zonal flow growth rate is rather weak. Hence, 
additional effects of negative magnetic shear, not included in the present 



Figure 4: The ZF growth (normalized to the TEM growth) as a function of 
magnetic shear s with normalized temperature gradient as a parameter. The 
results are shown for r\ e = 6.0 (squares) and r\ e = 4.0 (asterisk). The other 
parameters are k x = k y = 0.3, q x = 0.3, e = 0, rji = 0, e n = 1.0 and fraction 
of trapped electrons f t = 0.5. 
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study, may be needed in order to stabilize the turbulence. This observation 
may have some relevance for the observed difference in the dynamics of elec- 
tron and ion transport barriers in tokamak plasmas. More work, including a 
treatment of the coupled ITG/TE mode system using more realistic models, 
is needed in order to confirm this. This is left for future work. 
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